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Abstract 

We study the sensitivity of a 500 GeV e"'"e~ Linear Collider to CP 
violating anomalous couplings. We find that with 50 fb~^, and taking 
only one non-zero coupling at a time, the process e"'"e~ — > W~^W~ can 
be used to place the 95% confidence level bounds \R^\ < 0.1, \kz\ < 
0.1 and \gf \ < 0.1 from CP even observables. By studying certain 
distributions in the process e'^e~ — > ^~uv one of the bounds can 
be improved to \gf\ < 0.06. This process also allows the construction 
of a CP odd observable which can be used to place bounds on CP 
violating new physics. At the 95% confidence level we find Ik^I < 0.3, 
Ik^I < 0.2 and a much weaker bound for \gf\. 



1 Introduction 



In spite of its remarkable phenomenological success there are several aspects 
of the standard model that remain unexplained. Two of them are the mech- 
anism of electroweak symmetry breaking and the origin of CP violation. 
A broad class of models in which the electroweak symmetry is broken dy- 
namically by new strong interactions does not contain any new particles 
sufficiently light to be produced at a 500 GeV e~^e~ collider. New particles 
in these models have masses in the TeV range and only manifest themselves 
indirectly in experiments at lower energies. In general these models may 
violate CP and this would also manifest itself indirectly at low energy. 

It is convenient to describe the phenomenology of the most important 
features of this type of new physics in a model independent way. This is ac- 
complished by studying a low energy (below a few TeV) effective Lagrangian 
that contains only the standard model fields and where the effect of the new 
physics appears as higher dimension operators. These higher dimension op- 
erators modify the couplings of the observed particles, inducing "anomalous 
couplings" whose phenomenology has been studied in detail 

In this paper we study the effect of the lowest dimension operators that 
violate CP in the gauge-boson self-couplings. In has become standard to 
parameterize the three gauge boson coupling WWV following the notation 
of Ref. i: 

Cwwv = 9wwv{^-9lW+W-{d^V'' + d^V^ + -^^v^y.vo.pW^W^V'^f'^ (1) 

In writing this equation we have already dropped terms proportional to \y 
because they are of higher dimension. The overall normalization is gwW'^ = 
— e and gwwz = — ecot6'iy. Electro-magnetic gauge invariance forbids the 
term gj so we are left with three new CP violating parameters. One of them, 
g^ , violates C and conserves P; whereas the other two, and kz violate P 
and conserve C. 

The next to leading order electroweak chiral Lagrangian contains three 
CP violating operators whose couplings correspond to the lowest dimension 
contributions to the parameters in Eq. (|1|). They are |Q: 

^CP —odd 
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+ \a^igh>'^P''i:i{TW,,)i:i{TWp„) (2) 

o 

where we have used the notation of Ref . [Q . The relation between these cou- 
phngs and those in Eq. (|) is kz = e^{ai'i/ Cy^-aiA/ Sy^), = -(e^/s^)(ai3+ 

au) and = -eV(svFCw^)"i2 §• 

We do not wish to reproduce here all the details of the notation of Ref. [^], 
it is sufficient to remind the reader that the factor T that appears in all three 
operators constitutes an explicit breaking of custodial symmetry. This im- 
plies that if the new physics violates CP maximally (i.e. without suppressions 
from dimensionless parameters such as mixing angles) , the natural size of the 
coefficients is that of (f^/A^)Ap. These are the same arguments that give 
the natural size of the lowest dimension parity violating coupling 0. With 
V ^ 246 GeV and the scale of new physics A a few TeV we thus expect 
<^i2,i3,i4 ~ 10~^ if the symmetry breaking sector has a custodial symmetry to 
explain the smallness of Ap. In Ref. 0] Appelquist and Wu discuss a specific 
model in which they estimate that the coefficients ai2,i3,i4 are indeed of order 
10"'^ and correlated with Ap. On the other hand, if Ap is small accidentally, 
naive power counting tells us that these couplings could be at the few percent 
level, 0{v'^/A'^). These numbers will help us calibrate the significance of the 
constraints that we discuss. 

The best indirect bound that exists on any of these couplings is \k^\ < 
2 X 10~^, which arises from the neutron e.d.m. 0. This is a very tight 
constraint, but it is subject to naturalness assumptions. As usual, it is not 
a substitute for a direct constraint. Previous studies of W'j production at 
an upgraded Tevatron have concluded that it will be possible to place the 
constraint \k^\ < 0.1 [0]. This is of the same order as the bound that we find 
in this study for the 500 GeV e~^e~ collider, but a precise comparison is not 
possible without further knowledge of the experimental setups. 

There have been previous studies of the CP violating anomalous couplings 
in the process e~^e~ W~^W~ [§, but a detailed numerical analysis of 
the bounds that one can get at a NLC has not been done. The process 
e~^e~ — > uVZ has also been recently considered 0. The authors of Ref. 
find that one could place the bound g^ < 0.1 by studying a forward-backward 
asymmetry with 50 fb~^. 
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2 Bounds from observables in e^e W^W 

We start with the process e'^e" —>■ W^W~ at a center of mass energy of 
500 GeV without considering any specific decay channel for the W bosons. At 
this stage we also ignore the CP violating nature of the couplings k^, k^, gf 
and look only at the quadratic effects that they induce in the decay distribu- 
tion. It is possible to study truly CP- violating effects that are linear in the 
couplings in one of two ways. We could include absorptive phases in the form 
factors associated with Eq. (Q) 0, ^. If these phases arise from the same 
sector responsible for the anomalous couplings then they do not introduce 
additional suppression factors. This can be seen, for example, in the model 
of Ref. 0]. Alternatively, we could construct a CP odd observable involving 
the polarization vectors of the W bosons. This is equivalent to studying cor- 
relations that involve the momenta of the decay products of the W^-bosons 
to a specific channel. We take the second approach in the following section. 

In this section we consider the W bosons to be final state particles and 
take into account the efficiency for W^W~ pair reconstruction ey/w = 0.15 
[ p!0|] in our numerical simulation. Because we are ignoring, for now, the CP 
violating nature of the couplings, it is possible to bound them using the same 
CP even observables that we studied in Ref. The only difference be- 

tween the CP violating couplings that we study here, and the CP conserving 



couplings that we studied in Ref. |11| is that the CP violating couplings al- 



ways appear quadratically in the CP even distributions used to place the 
bounds.^ We use the same assumptions about systematic uncertainties and 
the same analysis of the differential distribution that we described in detail 
in Ref. In particular, we use a systematic error ~1.5%. This number 

arises from an uncertainty in the luminosity measurement of ~0.5%, an error 
in the acceptance an error for background subtraction ~0.5% and a 

systematic error in the knowledge of the branching ratio ~0.5%. From a 
analysis of da/dcos9 with 5 bins, we find the 95% confidence level bounds 
(taking only one non-zero coupling at a time): 

\k^\ < 0.1, \kz\ < 0.1, l^f I < 0.09 (3) 

The best bounds for this process are obtained using four bins, as discussed in 

^Again, this is because we are not including any possible absorptive phases. We can 
justify this a posteriori because the bounds that can be obtained from terms linear in the 
couplings and terms quadratic in the couplings are very similar due to the low statistics. 
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Ref. |TT[. However, the bounds we obtain using five bins are indistinguishable 
from those for four bins. We prefer five bins because this will be the optimal 
number for e"'"e~ — jj,^ jjL^vV and using the same number of bins in both 
cases will facilitate a comparison. 




Figure 1: Allowed (95% C. L.) region from the CP even angular distribution 
of e+e~ W'^W' for: a) kz — with = 0; b) k^ — with kz = 0; 
c) k^ — kz with gf = 0. 

In Figure |I| we present the allowed 95% confidence level regions when we 
take one of the CP violating couplings to be zero. We see that the bounds are 
indeed similar to those that can be placed at an upgraded Tevatron |0, and 
of the same order as the bounds that we found for CP conserving anomalous 
couplings |TT|. The fact that we obtain similar bounds for couplings that 
contribute linearly to the differential cross-section and for couplings that only 
contribute quadratically, already indicates that this process is not sensitive 
to the very small values predicted by naive dimensional analysis. 
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3 Bounds from observables in e^e fi^fi 

We now wish to consider a specific channel for the decay of the W^-bosons 
so that we can construct correlations that could single out CP- violating 
interactions. With this in mind, we need a final state that is easy to identify 
and that transforms into itself under CP. We thus choose to identify the W 
pairs by their /ii/ leptonic decays. We calculate the amplitudes for the process 
e+e^ fi^fi^uu and generate events for the three following subprocesses: 

Subprocess I — e^e^ — > ji^ ji^u^u^ , 
Subprocess II — e^e'^ ji^ ji^Ve^e , 
Subprocess III — e^e'^ fi^^^Ur^r ■ 

These subprocesses include 20, 21 and 11 Feynman diagrams respectively. 
The last one does not contain anomalous vertices, and constitutes pure back- 
ground. Sufficient events are generated with our Monte-Carlo to achieve a 
1% statistical error in the value of the cross-section. 

We first generate events with a cut on the muon scattering angle and on 
the muon pair invariant mass: 

170° <e< 10°, M^^ > 30 GeV (4) 

The angle 6 is the scattering angle between the fi~ and the e~ momenta in 
the e~^e~ center of mass frame. To study CP odd observables we need to 
make sure that our cuts are "CP blind", so the same cut is imposed on the 
angle between the and the e"*" momenta. These cuts are similar to the 
ones used by the LEP experiments. After imposing these cuts we assume a 
muon reconstruction efficiency equal to 1. 

The total cross-section for the e~^e~ fi^fi'up process with these cuts 
is 7.65 fb, which results (with an integrated luminosity of 50 fb~^) in only 
382 events. The resulting bounds will be limited by the small statistics so 
we will have to relax these cuts later on. 

^The current experiments at LEP have a typical region for muon reconstruction of 
170° < 9 < 10°. We assume that the experiments at a NLC will have a similar geometry 
so we use the same cut. The cut on the invariant mass of /x^/z" pair, m^^ > 30 GeV 
serves to reject Dalitz conversion of soft photons and to insure good angular separation of 
the muons. 
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We start by placing bounds on the CP violating anomalous couplings 
using the following observables: 

da da da 
dcosO ' dpf^ ' (ip^x/. ' 

where 9 is again the scattering angle between the /i~ and the e~ momenta 
in the e^e~ center of mass frame; is the muon three-momentum in the 
same frame, and p^x^i is proportional to the T-odd correlation ■ x p^-]. 
Numerically we work in the e^e~ center of mass frame and use: 

4 

Pmxm = ^Pt^-]z, (6) 

where the index z denotes the component along the beam direction. With 
this normalization p^x^l can take values from —1 to 1. If the polarization of 
the final leptons is not observed, and the beams are not polarized, this cor- 
relation serves to analyze the CP properties of the interaction |T2[. CP even 
interactions give rise to symmetric distributions (symmetric about the point 
P/ix/j = 0) in p^x^i, whereas CP odd interactions give rise to antisymmetric 
distributions in p^x/^- The antisymmetric distributions in this correlation 
will arise from interference between standard model amplitudes and the new 
CP-violating physics and will be linear in the new couplings. Notice that 
the correlation p^x^ is also odd under parity. This means that we will get 
terms proportional to k^^z from interference with the parity even standard 
model amplitude and a term proportional to gf from interference with the 
parity odd standard model amplitude. 

In Figure ^ we show the differential cross-section predicted by the stan- 
dard model at lowest order, as a function of cos6' and p^, for ^/s = 500 GeV. 
From these figures we see that the events predicted by the standard model 
are concentrated at small scattering angles and low muon momentum. Sim- 
ilarly we find that the standard model events have a symmetric distribution 
in p^x^t (as corresponds to CP conservation) that is very strongly peaked at 

Pf^Xf, = 0. 

In order to understand the effect of the cuts that we impose on the dis- 
tribution with respect to p^xfi it is convenient to express this correlation 
analytically. In addition to the angle 9 defined above, we need to define 9^, 
the angle between the yLi+ and the n~ momenta in the e^e~ center of mass 
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Figure 2: Standard model differential cross-section for the process e'^e~ — 
H~^fi~i'iy at y/s = 500 GeV with the cuts of Eq. (^) as a function of a) cos^ 
and b) p^. 



frame. Using a coordinate system with the z axis pointing in the direction of 
the /i~ momentum, the e~ momentum being in the x — z plane, and with 0^ 
being the azimuthal angle of the /i"*" in this coordinate system, the correlation 
is proportional to: 

P/,XM ~ sin sin 9f, sin 0^ (7) 

Working in the limit of massless muons, the angle 6fj_ is related to the invariant 
mass of the muon pair in the following way: 

We can now understand the effect of the cuts of Eq. (^: a) they remove 
the region of small sin 6' where p^xfi is small, thus increasing the signal to 
background ratio; b) they remove the region of small muon pair invariant 
mass effectively enhancing the region where cos 6^ is negative. 

In Figure ^ we show the deviations induced by the CP- violating couplings 
in the differential cross-section. For illustration purposes we use the values 

= 0.1, kz = 0.5, and = 0.5, with only one of them being non-zero 
at a time. The curves labeled '1', '2', and '3' correspond to the case of 
non-zero , kz, and respectively. It is clear from these figures that the 
kinematic regions where the new effects would be most important are high 
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muon momentum and backward scattering. The distribution with respect 
to p^x^t) is approximately symmetric about p^x^ = indicating that for the 
cuts in Eq. (^), the terms quadratic in the new couphngs (and thus CP-even) 
dominate. 




-0.4 -0.3 -0.2 -0.1 0.0 0.1 0.2 0.3 0.4 

Figure 3: Deviations in differential cross-section with the cuts of Eq. from 
their standard model values as a function of a) cos^^, b) and c) p^x^j.- In 
all cases the curves labeled 1, 2 and 3 correspond to gf = 0.1, kz = 0.5 and 
= 0.5 respectively. 

It is interesting to notice that for the set of cuts that we have used so far, 
the total cross-section is more sensitive to the value of than any of the 
distributions, we show this in Figure § 

To place bounds on the anomalous couplings, we used a standard x^- 
criterion to analyze the events (including the 0.5% anticipated systematic 
error in the luminosity measurement). We also investigated the sensitivity 
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Figure 4: Total cross-section with the cuts of Eq. (H) as a function of 
(sohd hne) and kz^-y (dashed-hnes indistinguishable in this figure). 

of the resulting bounds to different kinematic cuts and binning, but we did 
not find any way to enhance the sensitivity to the new couplings. This is 
probably due to the very low statistics available (382 events). We find that 
the best bounds are obtained by dividing the events into 5 bins. 

The results for y/s = 500 GeV, integrated luminosity of 50 fb~\ 5 bins, 
and at the 95% CL. are shown in Figure ^. In this figure the solid contours 
correspond to the bounds coming from the muon momentum distribution, 
the short-dashed contours from the scattering angle distribution, and the 
long-dashed contours from the correlation p^x^t- For k^^^z the best bounds 
arise from the muon momentum distribution, at about the same level as the 
bounds that this process places on CP conserving anomalous couplings 
For the bound is g^ < 0.06, slightly better than what we got in Eq. (|^), 
whereas for kz^-y the bounds are worse than those in Eq. @. 

Since our bounds are probably limited by the low statistics, we now study 
the effect of relaxing the cuts, and impose only the minimal cut: 

M^+^- > 5 GeV. (9) 

This is a very optimistic cut that still permits high efficiency in muon de- 
tection. Relaxing the cuts in this way has the effect of increasing the cross- 
section for the process e+e~ ji'^n^vu to about 113 fb. With an integrated 
luminosity of 50 fb~^ this results in 5660 events and consequently much better 
statistics. 
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Figure 5: 95% CL. bounds taking one coupling to be zero at a time and 
using the cuts of Eq. (Q). The sohd contours correspond to the bounds 
coming from the distribution, the short-dashed contours from the cos^ 
distribution, and the long-dashed contours from the p^xfj. distribution. 



In Figure |] we present the total cross-section as a function of the anoma- 
lous coupling (solid hne), kz (short-dashed line) and (long-dashed 
line). Comparison with Figure | shows that the relaxed cuts increase the 
sensitivity to the anomalous couplings. In Figure ^ we show the differential 
distributions with respect to the muon scattering angle, cos 9, and the muon 
momentum, p^. Once again we see that the standard model populates the 
regions of small scattering angle and low muon momentum preferentially. 
In Figure |^ we show the change in the differential cross-section when the 
anomalous couplings take values = 0.1, kz = 0.5, and = 0.5. We 
take only one non-zero anomalous coupling at a time and use these values 
for illustration purposes only. The standard model distribution is largest 
near 6 = whereas the new physics contributions are largest near 6 = 90°. 
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Figure 6: Total cross-section with the cuts of Eq. as a function of gf 
(sohd hne) and ^^,7 (dashed-hnes indistinguishable in this figure). 

Nevertheless, we find better sensitivity to the new physics when we do not 
impose the angular cut that excludes the region of small 6*, indicating that 
our analysis is limited by statistics. 

To place bounds on the anomalous couplings, we use a standard x^- 
criterion to analyze the events, include the anticipated 0.5% systematic error 
in the luminosity measurement, and take the muon identification efficiency 
to be 1. We find that the best bounds are achieved by dividing the events 
into 5 bins. The 95% C.L. results for ^/s = 500 GeV, integrated luminos- 
ity of 50 fb~^ and 5 bins for the bounds on gf , kz, and shown in 
Figure |^. In this figure we have taken one of the three couplings to be zero 
and looked at the projection of the allowed region into the plane of the other 
two couplings. In this figure the solid contours correspond to the bounds 
coming from the correlation p^x/^, the short-dashed contours from the muon 
momentum distribution, and the long-dashed contours from the scattering 
angle distribution. Comparing Figures ^ and ^ one can see that with the 
relaxed cuts, the best bounds are obtained from the correlation p^x/j,- The 
improvement in the bounds is partly due to the increased statistics, and 
mostly due to the fact that when we relax the cut on the muon pair invariant 
mass we include a region of phase space that has the largest sensitivity to 
Pf_ix^i- In view of Eqs. (|^) and (^), the region of smaller muon pair invariant 
mass (M^+^- < 30 GeV) appears to contain the region where sin is large. 

Alternatively, taking only one non-zero coupling at a time we find (from 
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Figure 7: Standard model differential cross-section for the process e'^e~ 
jji'^H'i'V at a/s = 500 GeV and the cuts in Eq. (||) as a function of a) cos^ 
and b) p^. 



the correlation of Eq. 



\k^\ < 0.27, \kz\ < 0.18, \gi\ < 0.0^ 



(10) 



4 Bounds from a CP-odd observable 

In the previous section we have seen that the bounds that can be placed on 
the anomalous couplings using the correlation p^x^ are of the same order as 
those that can be placed from other observables. Therefore, it is interesting 
to see whether one can isolate the CP-odd components of the distributions 
with respect to p^x^, and in that way be able to really bound new CP 
violating interactions. 

To understand the effect of the different sets of cuts on these bounds we 



present in Figures |T0| and |TT| the differential cross-section with respect to the 
correlation p^x/x for ^-y = 0.5 and kz = 0.5 respectively. 

In these figures we have separated the contributions to the differential 
cross-section arising from terms linear in the anomalous couplings (the curves 
that are antisymmetric about p^xfi = 0) from those arising from terms 
quadratic in the anomalous couplings (the curves that are symmetric about 
P/ixfi = 0). We also show how these results vary when we go from the stronger 
cuts of Eq. (H) to the relaxed cut of Eq. (0). 
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Figure 8: Deviations in differential cross-section with the cuts of Eq. from 
their standard model values as a function of a) cos 9, b) and c) the CP odd 
correlation of Eq. (H). In all cases the curves labeled 1, 2 and 3 correspond 
to (?f = 0.1, = 0.5 and = 0.5 respectively. 



Notice that the normalization of the curves with strong and weak cuts is 
different as it corresponds to the respective total cross-section. From these 
curves we see that the relaxed cuts are not only better because they increase 
the statistics, but they also increase the relative contribution of the truly 
CP-odd term linear in k^^z as we argued in the previous section. 

In order to quantify the bounds that can be placed in this way, we intro- 
duce the integrated CP-odd observable 

^ = / ■ sign(p^xM)%xM- (11) 
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Figure 9: 95% CL. bounds taking one coupling to be zero at a time and 
using the cuts of Eq. (|). The short-dashed contours correspond to the 
bounds coming from the distribution, the long-dashed contours from the 
cos 6' distribution, and the solid contours the P/j^xn distribution. 



Specifically we take the sign(p^xAt) to be zero if p^xfi = to exclude that 
point, and use the following criterion to place bounds; 

I / 7—— ■ sign(p^xM)'^PMXMl < 2 ■ Aaexp ■ (12) 

The right-hand side of this equation corresponds to two standard deviations 
for the experimentally measured cross- sect ion. Given the definition of A it 
is appropriate to use the same experimental uncertainty of the total cross- 
section: 

syst + Slat , (13) 
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Figure 10: Deviations in the differential cross-section with a) the cuts of 
Eq. and b) the cuts of Eq. (|^) as a function of p^x^i 

for = 0.5. We 

have separated the contribution from the term hnear in (antisymmetric 
curve) from that due to the term quadratic in k^ (symmetric curve). 



where ^ 

8 Stat = I ■ (14) 

The best bounds come from using the relaxed cut of Eq. and are given 

\~kz + 0.7 ■ ~k^\ < 0.2 (15) 
It is also possible to obtain bounds on but they are much weaker. 



5 Conclusions 

We have studied the effect of CP violating anomalous couplings on the pro- 
cess e~^e~ —>■ W^W~ . Using CP even observables we have found that a 
NLC with ^/s = 500 GeV and 50 pb~^ can place the bounds \k^,z\ < 0.1 and 
\gi\ < 0.09. These bounds are comparable to those that can be placed with 
an upgraded Tevatron, and are of the same order as the bounds that can be 
placed on CP conserving anomalous couplings. These bounds originate in 
the quadratic contributions of the couplings to the differential cross-section. 
By looking at the fiu decays of the ly-bosons we were able to construct a 
CP odd correlation that can directly bound the CP-violating terms (linear 
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Figure 11: Deviations in the differential cross-section with a) the cuts of 
Eq. (H) and b) the cuts of Eq. as a function of p^x^t for kz = 0.5. We 
have separated the contribution from the term hnear in kz (antisymmetric 
curve) from that due to the term quadratic in kz (symmetric curve). 



in the couphngs) in the differential cross-section. We found that it will be 
possible to place the bounds \k^\ < 0.3 and \kz\ < 0.2. We conclude that 
the sensitivity of a NLC to CP violating anomalous couplings is similar to 
its sensitivity to CP conserving anomalous couplings. From our dimensional 
analysis we also conclude that it is unlikely that CP violating anomalous 
couplings will be seen by a NLC unless the smallness of Ap is accidental. 
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